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Abstract 

We perform a detailed analysis of the one-loop corrections to the light neutrino mass matrix 
within low scale type I seesaw extensions of the Standard Model and their implications in ex¬ 
perimental searches for neutrinoless double beta decay. We show that a sizable contribution to 
the effective Majorana neutrino mass from the exchange of heavy Majorana neutrinos is always 
possible, provided one requires a fine-tuned cancellation between the tree-level and one-loop con¬ 
tribution to the light neutrino masses. We quantify the level of fine-tuning as a function of the 
seesaw parameters and introduce a generalisation of the Casas-Ibarra parametrization of the neu¬ 
trino Yukawa matrix, which easily allows to include the one-loop corrections to the light neutrino 
masses. 
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I. INTRODUCTION 


The main goal of this work is to study in detail under which conditions the right-handed 
(RH) neutrinos present in a general type I seesaw scenario [T] can give a direct sizable 
contribution to the neutrinoless double beta (Ozz/3/3) decay rate, i.e., a contribution in the 
range of sensitivity of the current and upcoming 0 v(3/3 decay experiments, once all the 
relevant constraints are included in the analysis. 

In [21S], it was shown that a sizable sterile neutrino contribution to the 0 vf3f3 decay can be 
achieved if the heavy neutrino spectrum is hierarchical, with at least one RH neutrino with 
mass M below 100 MeV and the other state(s) above this scale. However, this spectrum is 
disfavoured by cosmological observations since the region M e [1 eV, 100 MeV] is excluded 
by BBN and CMB data [H Ej. In [6H8] the possibility of having a relevant contribution 
from heavy RH neutrinos up to the TeV scale was explored. 1 It was found that indeed RH 
neutrinos as heavy as 100 GeV-10 TeV could, in principle, give a sizable and observable 
contribution to the 0 vf3j3 decay rate. In [S] the role of the fine-tuning and one-loop effects 
were discussed, concluding that for RH neutrino masses above 10 GeV a relatively high level 
of fine-tuning would be required. In [3j a more detailed study of the one-loop effects was 
performed and it was found that indeed they are significant and can play a very important 
role in the type I seesaw scenario. The lepton number violation introduced through the RH 
neutrino Majorana mass term, required to obtain a sizable effect in the 0zz/3/3 decay rate, 
naturally appears at one-loop level in the light neutrino sector. If fine-tuning is not invoked, 
the light neutrino mass constraints on the one-loop corrections make it very difficult to 
obtain a significant (RH) heavy Majorana neutrino contribution in the Cffi/3/3 decay effective 
Majorana mass, i.e., to have |m^ avy | > 0.01 eV, mj^ avy being the heavy Majorana neutrino 
contribution under discussion. We will show, in particular, that the scenario in which RH 
neutrinos with a mass M > 1 GeV can give a sizable contribution to the 0 u/3/3 decay 
rate necessarily involves a fine-tuned cancellation between the tree-level and one-loop light 
neutrino contributions. 

More specifically, in this work we re-analyse the conditions under which the heavy Majo¬ 
rana neutrinos with masses M > 100 MeV of the type I seesaw scenario can give a significant 

1 The interplay between the light and heavy Majorana neutrino contributions in 01 // 3/3 decay was investigated 
phenomenologically first in [9]. 
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direct contribution to the Ou/3/3 decay effective Majorana mass, i.e., a contribution in the 
range of sensitivity of the current and upcoming 0 u/3/3 decay experiments. We show that 
for M > a few GeV this requires a relatively large active-sterile neutrino mixing (charged 
current couplings of the heavy Majorana neutrinos). We clarify which seesaw realisations 
can provide the requisite mixing. We discuss the impact of the one-loop corrections in the 
different type I seesaw realisations considered. We analyse also numerically the problem 
of the sizable heavy Majorana neutrino contribution to the 0 u/3/3 decay effective Majorana 
mass, by studying the full parameter space, including the relevant one-loop corrections and 
the bounds on the active-sterile neutrino mixing from direct searches, charged lepton flavour 
violation and non-unitarity PUHZ!. We quantify, in particular, the level of fine-tuning re¬ 
quired in order to have a sizable heavy neutrino contribution to the 0 u/3/3 decay rate. In 
order to do the analysis and generate the right pattern for the light neutrino masses and 
mixing, we have constructed a modification of the Casas-Ibarra parametrization [18] , which 
takes into account the impact of the one-loop corrections. 

The paper is organized as follows: in section [TT] we derive under which conditions it is pos¬ 
sible to obtain a sizable active-sterile neutrino mixing, which can strongly affect the effective 


Majorana neutrino mass, mpp. In section III we study the impact on rnpp of the one-loop 
corrections to the light neutrino masses and present our modified Casas-Ibarra parametriza¬ 


tion which takes into account the one-loop effects. In section IV we perform the numerical 
analysis and quantify the level of fine-tuning necessary to have a dominant contribution in 
rriQp from the exchange of the heavy (sterile) neutrinos. Finally, we summarise our results 
in the concluding section. 


II. LARGE ACTIVE-STERILE NEUTRINO MIXING AND Ou/3/3 DECAY 

We consider the most general type I seesaw scenario [1] with n > 2 RH neutrino fields 
u sR (s — 1, ■ ■ ■ ,n). After the spontaneous breaking of the electroweak (EW) symmetry the 
full neutrino mass Lagrangian is 


C v = - u tL (rn D )i s u sR - - u c sL (M R ) st u tR + h.c. 


( 1 ) 


where (. = e, /r, r and Kl = C u sR T , C being the charge conjugation matrix. M R = ( M R ) T 
is the Majorana mass matrix of the RH neutrinos and m R is the 3xn neutrino Dirac mass 
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matrix. The full mass matrix derived from Lagrangian (|Tj) is therefore 


M 


( O m D 

m T D M r 


U* diag (rrii, M k ) U\ 


( 2 ) 


where rrii (i = 1,2,3) and M k (k = 1 , n) are the light and heavy Majorana neutrino 
masses, respectively. We define O as a 3 x 3 matrix with all elements equal to zero. The 
full neutrino mass A4 is diagonalised by a (3 + n) x (3 + n) unitary matrix U, through a 
well known rotation between the neutrino flavour and mass eigenstates. We give below the 
relation between the left-handed (LH) components of the corresponding fields {yn,, v c s l an( i 
XiL, N kL )\ 


VlL 


= U 


V, 


sL , 



(3) 


Taking into account that the active block of U is unitary to a very good approximation, the 
complete mixing matrix can be expanded as 2 


U 


( 1 - 00t/2 

V 


1 - 0 + 0/2 


U PMNS 0 j + 0 (fl3) 


( h^PMNS @V 
— 0lf/pMNS V 


+ ^ (0 2 ) , (4) 


where 6 is a 3 x n matrix with “small” entries, which characterises the mixing between the 
active and the sterile neutrinos, ITpmns is th e PMNS neutrino mixing matrix mm and V 
is a n x n unitary matrix. The quantity (6 V)i k , £ — e, /r, r, k — 1,..., n, is the coupling of 
the heavy Majorana neutrino N k to the charged lepton P in the weak charged lepton current, 
and to the flavour neutrino vi in the weak neutral lepton current. 

From the diagonalization of the complete neutrino mass matrix A4, at leading order in 6 
we have [B] 


6 M r 0i ~ ^ PMNS rh f/p MNS , 
6* M r « rn D , 

M r » V* M V f , 


(5) 

( 6 ) 
(7) 


where 


rh = diag(mi, m 2 , m 3 ), M = diag (Mi,..., M n ). 


( 8 ) 


2 In the following we work in the basis in which the charged lepton mass matrix is diagonal. 
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It follows from Eqs. ([5]) and (|7j) that 


(evy m (evy ^ -u; MNS muj 


PMNS 


(9) 


In terms of the seesaw parameters we have for the active-sterile neutrino mixing: 


9* ~ mo M r 


-l 


( 10 ) 


Using Eqs. ([5]) and (10), we recover the usual type I seesaw relation for the (tree-level) light 
neutrino mass matrix, namely 


= -rn D M R ] rn T n = -9* M R 9 ] = —(9 V)* M (9 V)^ = E/£ MNS m Up 


PMNS ' 


( 11 ) 


The effective Majorana neutrino mass, mpp, which enters in the 0 u/3/3 decay amplitude, 
receives, in general, two different contributions, corresponding to the exchanges of the light 
and heavy virtual Majorana neutrinos: 


light - heavy 

m 99 = + m pp , 


with 


m 


light _ 

99 ~ 


(flpMNs)e 




( 12 ) 


(13) 


i=l 


where we have used Eq. (J9J), which holds at tree-level in the type I seesaw models. A good 
estimate for the contribution due to the heavy Majorana neutrino exchange for Mk 100 
MeV is 0 

(14) 


rn 


heavy 

99 


(MJMtfMt, 


where M a ~ 0.9 GeV and f(A ) depends on the decaying isotope considered. For, e.g., 48 Ca, 
' 6 Ge, 82 Se, 130 Te and 136 Xe, the function f(A ) takes the values f(A) ~ 0.033, 0.079, 0.073, 
0.085 and 0.068, respectively. 


Using Eq. (14), it is easy to estimate the minimum mixing (#U) m in required in order to 
have a contribution at the aimed sensitivity of the next generation of 0 u/3/3 decay experi¬ 
ments, that is |m^ avy | > 10 -2 eV. In Fig. 0 we compare this estimate for (^U)^ for the 
' 6 Ge isotope, (0U) 2 lin ~ 1.6 x 10~ 10 M GeV^ 1 (dashed line), with the naive seesaw scaling 


suggested by Eq. (11), (#U) 2 aive = \J Am atm /M ~ 5x 10 11 GeV /M (solid line) as a function 
of the RH neutrino mass scale M (expressed in units of GeV). 


5 







FIG. 1: Active-sterile neutrino mixing. The dashed line stands for an estimate of the 
minimum ( 6V ) 2 required in order to have |m^ avy | > 10~ 2 eV in the case of Ou/3/3 decay of ,6 Ge. 
The solid line corresponds to the naive seesaw scaling of (9V) 2 (see the text for further details). 

From Fig. [l] it is clear that for RH neutrino masses larger than ~ 1 GeV a considerable 
enhancement with respect to the naive seesaw scaling of 6V is required in order to have a 
sizable RH neutrino contribution. Obviously, this enhancement increases with the mass of 
the RH neutrinos. We notice that in the region M ps 500 MeV-1 GeV, the naively estimated 
mixing, (6W) 2 aive , is in the right ballpark. Similar conclusions are valid for (6W) 2 nin and 
(6W) 2 aive in the cases of Ou/3/3 decay of other isotopes ( 48 Ca, 82 Se, 130 Te, 136 Xe, etc.). 

A. Casas-Ibarra Parametrization and Large Active-Sterile Neutrino Mixing 

In order to understand under which conditions an enhancement with respect to the 
naive scaling of the active-sterile mixing (or equivalently, of the charged current couplings 
of the heavy Majorana neutrinos ( 0V)ik ) can be expected, we employ the Casas-Ibarra 
parametrization of 9V DSI- In this parametrization the light neutrino masses and the angles 
and phases of the PMNS matrix are input parameters, in such a way that the correct light 
neutrino mixing pattern is always recovered. The Casas-Ibarra parametrization is obtained 
rewriting Eq. ([5]) as 

(±i m " 1/2 4 mns W M 1/2 ) (±im" 1/2 4 mns 6V M 1/2 ) T = RR T = 1, (15) 
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where R is a general 3xn complex matrix which parametrizes the new physics degrees of 
freedom associated to the sterile neutrino sector. Using this parametrization, QV can be 
written as 

QV = T* U PM Nsm 1/2 RM~ 1/2 . (16) 

The matrix V can be set to the unit matrix if one works in the basis in which the Majorana 
sub-matrix Mr is diagonal. 3 


Naively, from Eq. (11) one may conclude that QV ~ O ( he., that the mixing (or 


coupling) QV is expected to be suppressed by the heavy neutrino mass scale. However, having 
a larger mixing is perfectly possible due to an enhancement factor contained in the matrix 
R 0IZ|. Obviously, such enhancement can only be in agreement with the light neutrino 
spectrum if there is a non-trivial suppression/cancellation in the l.li.s. of Eq. ([9j . This 
extra suppression is related to particular textures of the neutrino mass matrix, which can be 
motivated, for instance, introducing an extra [/(1) global symmetry in the Lagrangian, as it 
is the case in the so called “inverse” and “direct” seesaw models [231 123] , In these models 
the indicated global symmetry can be identified with that corresponding to the conservation 
of a non-standard lepton charge (see further). 

In the following we will focus on the minimal seesaw scenario with n = 2 RH sterile 
neutrinos 4 (see, e.g., ESI) giving rise to two heavy Majorana mass-eigenstate neutrinos, 
which predicts one massless and two massive light active neutrinos. For the light neutrino 
mass spectrum with normal hierarchy (NH) and inverted hierarchy (IH) we have 


m 1 = 0, m 2 = \j A m 2 2l , m 3 = sj A m 2 3l , (NH) 

mi = \J | Am| 2 | — Amjj, m 2 = \J\ Am| 2 |, m 3 = 0 

The current best fit values obtained from the global fit analysis in [2SJ are 


(IH). 


(17) 

(18) 


A m 2 2l = 7.50 x 10" 5 eV 2 , (19) 

A m 2 3l = 2.457 x 10“ 3 eV 2 (NH) and Am 2 2 = -2.449 x 10" 3 eV 2 (IH) . 


In this minimal seesaw scenario, the two (tree-level) non-zero light neutrino masses m 2 


3 An extension of this parametrization to all orders in the seesaw expansion can be found in 21, 22 . 

4 In the present article we will use the term “heavy Majorana neutrinos” for Majorana neutrinos having 
masses exceeding approximately 100 MeV. 
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and m| ree (m^ ree ) in the case of NH (IH) neutrino mass spectrum satisfy the relation: 


m 2 reem 3 W = — detfM^ 1 ] detf-mjm^i], NH (IH), 


( 20 ) 


which is basis independent. 


In the considered case the /?-matrix, which enters into Eq. (16), can be parametrized as (7] 

/ n \ 


R = 


0 0 

cos (0 45 + iy) - sin (0 45 + iy) 
^sin (0 45 + iy) cos ( 0 4 5 + i'y) ) 

(cos (0 45 + i'y) - sin (0 45 + iy)^ 


for NH, 


( 21 ) 


R = 


sin (0 45 + i'y) cos (0 45 + iy) 

0 0 


for IH, 


( 22 ) 


where 0 4 5 and y are real parameters. If R were real, i.e., y = 0, there is no way to obtain 
any enhancement of the couplings/mixings 6V of interest since R would essentially be a real 
orthogonal matrix. However, for y ^ 0 and e ±7 1 an enhancement of 6V is possible: 


| cos (0 45 + iy) |"’ = cos 2 0 45 + sinh 2 y > 1<S> e ±7 1, 
| sin (0 45 + iy) | 2 = sin 2 0 45 + sinh 2 y » 1 <=> e ±7 1. 

In fact, for e ±7 3> 1 the expression of R in the NH case reduces to 


(23) 


R « e~ id45 — 


3±7 


^0 0 ^ 

1 ±i 

V* l J 


NH. 


(24) 


Similarly, one can derive from (22) the same limit of R for the IH neutrino mass spectrum: 

/i _u-\ 


R « e“ i04S — 


=±7 


1 ±i 
Ti 1 
0 0 


IH. 


(25) 


Notice that the Casas-Ibarra parameter y in (24) and (25) can be related to the maximum 
eigenvalue y [Tj of the Dirac mass matrix rnu in Eq. (J2]), that is 

1 


y 2 v 2 = 2 max jeig {mo mj^ | = -e ±7 Mi (m 2 + m 3 ) (2 + z) , 
y 2 v 2 = 2 max |eig (mp m^ | = -e ±7 M 1 (mi + m 2 ) (2 + z) , 


NH 

IH. 


(26) 

(27) 














where z denotes the relative mass splitting of the two heavy Majorana neutrino masses, 
z = (M 2 — Mi)/Mi, and v = 246 GeV is the EW symmetry breaking scale. 


Introducing the expression (24) (or (25)) in Eq. (16) one obtains [51 171 1271 


m 


n 


, . M 2 

±m/ mT 


Then, in terms of y the active-sterile neutrino mixing in Eq. (16) takes the form [7] 

2 


n\ 


2 _ 1 y 2 v 2 m 3 

U e3 + iy/m 2 /m 3 U t 2 

2 (2 + z) Ml m 2 + rn 3 
|2 1 y 2 v 2 m 2 

2 (2 + z) Ml mi + m 2 

Ue 2 + *V m i/ m 2 Uei 


NH 

IH. 


(28) 

(29) 

(30) 


All in all, the previous relations imply that in the basis in which the RH neutrino Majorana 
mass term is diagonal, the neutrino Yukawa couplings, or equivalently (mo)n and (moje, 
should satisfy the following relation: 

(m D )ei 


, - Mi 

±S,/ iV4 


(31) 


(™ D ) i2 

Any texture of the neutrino mass matrix which satisfies this condition gives rise to relatively 
large couplings 6V with the right suppression/cancellation in the light (flavour) neutrino 
mass matrix, which allows to recover the correct light neutrino mass spectrum at tree-level. 
The relatively large QV thus generated can saturate the present bounds even in the case in 
which the heavy Majorana neutrino spectrum is hierarchical. 


Using Eqs. (14) and (28), one can easily estimate the contribution to the Ov/3/3 decay 
effective Majorana mass due to the exchange of the heavy Majorana neutrinos in the large 
coupling/mixing case of interest [6j: 

\ 2 ' 


heavy 

l /3/3 


ma'a'" ~ - (9V) el f(A) -f- 



1 - 


Ml 


Mi + AM 


(32) 


with 5 AM = M 2 — Mi. Clearly, if AM <C Mi the contribution will be proportional to AM, 
while in the limit AM Mi the dependence on AM is subleading since the lightest RH 
neutrino dominates the contribution. 

The interplay between the light and heavy Majorana neutrino exchange contributions in 
the effective Majorana mass, rnpp = + m^ avy , in the scheme under discussion in which 


5 Note that (0V) 2 el depends, in particular, on the phase $ 45 - This implies that mpp will also depend on 
045 [ 7 ;. 
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Eq. (28) holds and mj^ avy is given by Eq. (32), was investigated in detail in [7] in the case 
when the two heavy Majorana neutrinos form a pseudo-Dirac pair, 0 < AM = M 2 — M\ <C 
Mi, M 2 , and have masses in the interval ~ (50 — 1000) GeV. It was found that there exists 
a relatively large region of the allowed parameter space of the scheme in which the heavy 
Majorana neutrino contribution can change drastically the predictions based on the light 
Majorana neutrino exchange contribution. More specifically, it was found that [7]: i) \mpp\ 
in the case of NH spectrum can have values in the interval 0.01 eV < [m^) < 0.1 eV, i.e., 
in the range of sensitivity of the current GERDA [29], EXO [30], Kamland-Zen ]7l! and 
CUORE [32] experiments and of a few other experiments under preparation (Majorana [53], 
SNO+ [33], AMORE [35], etc.). We recall that in the case of Oz/d [3 decay generated only by 
light Majorana neutrino exchange we have (see, e.g., [36] [37]) [m^l = |mj]| ht | < 0.005 eV; 
ii) [m^l in the case of IH spectrum can be strongly suppressed due to partial, or even total, 
cancellation between and m^ avy in mpp (see also [5%]). Since the magnitude of m^ avy , 
as it follows from Eq. (32), depends on the atomic number A of the decaying nucleus [9], the 
cancellation between and m^ avy in rrigp can take place for a given nucleus (say, e.g., 

for 48 Ca) but will not hold for other nuclei ( 76 Ge, 82 Se, 130 Te, 136 Xe, etc.). If the 0 u/3/3 decay 
is due only to the light Majorana neutrino exchange we have in the case of IH spectrum, as 
is well known [39] (see also, e.g., [36]), 0.013 eV < [m^l = < 0.050 eV. 

On the other hand, in [7| the role of the one-loop corrections was not studied. In [3] it was 
shown that the one-loop corrections to the light neutrino masses generated in the scheme 
under discussion turn out to be very relevant. Essentially, a sizable heavy contribution to 
the 0 u/3/3 decay for heavy masses in the range ~ (50 — 1000) GeV generates at the same 
time a very large one-loop correction to the light neutrino masses. In this work we analyse 
in detail the role of the one-loop effects showing that similar conclusions to the ones drawn 
in [7] will be obtained. However, we will also show that the price one has to pay in order to 
have a significant impact of the heavy neutrinos in the 0 u/3/3 decay is the requirement of a 
highly fine-tuned cancellation between the tree-level and one-loop contributions to the light 
neutrino masses. 
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B. Comparison with Extended and Inverse Seesaw Scenarios 


As an application of the previous results, we consider the effect of heavy RH neutrinos 
on the 0 u/3/3 decay amplitude in the case of two different realisations of the type I seesaw 
scenario, which predict a large active-sterile neutrino mixing 6V, that is the well known 
extended seesaw (ESS) |3D] and inverse/direct seesaw (ISS) [23, 23] models. In particular, 
we will clarify how the large mixing realisations described in the previous section in terms 
of the Casas-Ibarra parametrization match with the ISS and ESS scenarios. 

In order to understand the predictions in these classes of models it is useful to adopt the 
following parametrization of the generic mass terms in the seesaw Lagrangian 0 - namely 


M = 



( 


O Ym/x/2 eY 2 v/V2 


\ 


Yjv/V2 

eY lv/y/2 


A 


A 

R 


(33) 




where Y\ = (y ie , y^, yi T ) T , for i — 1,2. This parametrization is completely general and, in 
principle, e, y, y! and A can take any value. 6 However, e, y and y' can be interpreted as lepton 
number violating couplings and, therefore, in principle they take arbitrarily small values, 
because in this case there is an approximate global symmetry of the seesaw Lagrangian 
corresponding to the conservation of the lepton charge L' = L e + + L r + L\ — L 2 , where 

L 1 and L 2 are the charges carried by the RH neutrino fields vm and z/ 2 r, respectively. In 
the limit of e = y = y! = 0, the conservation of LI is exact. In this case the neutrino sector 
consists of three massless neutrinos and one massive Dirac fermion, which can be inferred, 
in particular, directly from the expression of the charge LI in terms of the charges Lp and 


l i )2 [mug. The exact conservation of L' corresponds to the case in which condition (t28|) 
is exactly fulfilled and the RH neutrino splitting satisfies: AM = M 2 — Mi —> 0. 

In terms of the new parameters, the exact (tree-level) expression of the light neutrino 
mass matrix given in ([TT|) is proportional to y and e, that is 


m* ree = 


(y Yi Yf + e 2 y' Y 2 Y 2 T - A e (Y 2 Yj + Y, Y^)) , (34) 


2 (A 2 — y'y) 

and thus if y = e = 0 there is a complete cancellation at tree-level for the light neutrino 
masses. As we will see in the next section, if y! is different from zero, at least one neutrino 


In the following we will assume for simplicity that all the parameters introduced in Eq. (331 are real. 
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mass can be generated at one-loop, even for fi = e = 0 [3]. Furthermore, from the diago- 


nalization of (34), we obtain for the product of the smallest (rnf ee ) and the largest (m^ ree ) 
light neutrino masses: 


(35) 


1 tree tree 1 


\m l m h 

= | 


de 2L2 


{yin + Vlr) + yleivln + vlr) - Zy^eivi^ix + yiry2r) + (?/2 M 2/lr “ 2/l/d/2r) 2 | 

4| A 2 — nn'\ 

From this relation it follows that in order to have two massive active neutrinos at tree-level, 
i.e., m y 0, i ) an explicit breaking of the lepton charge conservation via the neutrino 
Yukawa couplings is necessary, that is the parameter e must always be different from zero; 
ii) the vectors of neutrino Yukawa couplings Yi and Y 2 cannot be proportional. 


Accordingly, the two seesaw limits of Eq. (34) which give rise to large active-sterile 


neutrino mixing 9V and generate sufficiently small active neutrino masses are: 

• i) /T A, y\ a v /U, ey 2a v (ESS limit). This limit matches the so-called extended 
seesaw [40J models and corresponds to a hierarchical spectrum for the heavy neutrinos: 

A / A 2 


M 1 « (A 2 /// - n ), (8V) 


n 


V2Mi [ VU n'-n 


M 2 « n' + A 2 ///, (8V) 


12 


y/2M 2 


yu i 


e V2t ( 1 


A 2 


2(/i'-/i) 2 


2(/i'-/i) 2 7J ’ 
(36) 
A 


+ ey 2 £ — 


H' — n 


(37) 


where we also show the corresponding mixing with the active neutrinos. Then, the 
approximate tree-level contribution to the 0 u/3/3 decay effective Majorana mass due to 
the exchange of the light and the heavy neutrinos is 

,2 ' - A 


m 


light 

73/3 


2 (A 2 //i' — n) \n' 2e fi' yieV2 ‘ 


m 


heavy 

w 




v 2 M 2 


/A 2 


A 


2 (A 2 /m'- A*) 3 U' 2 y ‘ e 2 %' aie !,2e 


(38) 

(39) 


respectively. The dominant term in m^ avy is due to the exchange of the lighter of the 
two heavy Majorana neutrinos N \, the exchange of N 2 giving a subleading (and negli¬ 
gible in the leading approximation we employed) correction. Notice that, if A 2 /fj /i, 
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m^| ht becomes independent of y! while m^ avy is proportional to y 


/ 3/3 


m 


light 

PP 


2 A 2 


(h A? e - 2 e A ?/i e r/ 2e ) , 


m 


heavy 


n! v 2 M 2 

i(A)L yM 


Vie ~ 2f jVleV2e 


(40) 

(41) 


ii) A |/i a n //, y, e y 2a v (ISS limit). This limit corresponds to a minimal realisa¬ 
tion with only two RH neutrinos of the so-called inverse or direct seesaw models [ 28) . 
In this case the heavy neutrino spectrum is quasi-degenerate, forming a quasi-Dirac 
pair 


M 1 


Mo 


A- 


A 


A + h' 


A + A 


m 


a 


(0V) 


a 


2Mi 


2Mo 


An 1 + 


A- A 
4A 


- e Z/2/ 1 - 


A-A 
4A 


Aw 


A — A 
4A 


+ eA« 1+ 


A — A 
4A 


(42) 


(43) 


In this limit the light and heavy contributions to the 0uj3/3 decay rate are given by: 


m 


light 

PP 


2 A 2 


(AA? e -2eAy le y 2e ) , 


m 


heavy 


V 2 M 2 

f( A )~ t/(( 2 A + a')aL - 2eAy le y 2e ) 


(44) 

(45) 


W " J v “ / 2 A 4 

Both of them are proportional to the small lepton number violating parameters, as it 
should be. Notice that the expression of m^| ht above is exactly the same as the one 


given in Eq. (40). 


On one hand, it follows from Eqs. (40), (41), (44) and (45) that a relatively large con¬ 
tribution to the 0 u/3/3 decay rate due to the heavy Majorana neutrino exchange might be 
possible at tree-level without affecting the smallness of the light neutrino masses since in 
the limits considered here m^ avy oc y\ while is independent of y'. On the other 

hand, Eqs. (36 37) and (42 43) confirm that the condition to obtain relatively large mix¬ 


ings, Eq. (28), is fulfilled at leading order, that is in the Casas-Ibarra parametrization the 


A-matrix corresponding to these two cases is similar to the textures reported in Eqs. (24) 


and (|25j). 

Finally, we note that in the case of the ISS model, the smallness of the light neutrino 
masses comes from the existence of an approximate symmetry corresponding to the conser¬ 
vation of the lepton charge L'. In contrast, in the ESS limit, the conservation of the lepton 
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charge L' is strongly violated through the [j! coupling. This means that, in principle, the 
one-loop corrections to the neutrino masses can be expected to be more important in the 
ESS limit than in the ISS one since in the ESS case there is no symmetry protecting the 
light neutrino masses from getting relatively large corrections [3]J. 


III. ONE-LOOP CORRECTIONS TO THE NEUTRINO MASS MATRIX 


We turn now to the computation of the one-loop corrections to the light neutrino mass 
matrix and the effective Majorana neutrino mass associated to 0 v(3/3 decay amplitude. 

At one-loop the neutrino self-energy E(p) provides the dominant finite correction to 
m v ||3j 05H18], which depends on the square of the neutrino Yukawa couplings, as in the 
tree-level contribution 0- and is further suppressed by the one-loop factor 1/(16 n 2 ). In 
a generic basis, with the Dirac and Majorana mass terms defined in Lagrangian (JT|) , we 
obtain: 


j 1—loop 

M = \ I =f/*diagK,M fc )f/ t , 

m T D M r 

where the new Majorana mass term generated at one-loop is in this case 

1 


(46) 


m 


1—loop _ 


(4 7T v) 2 

The loop function F(x) is defined as 


m D ( Mr 1 F(M r Mr) + F(MrMr) Mr 1 ) mj, 


t 


(47) 


F(x ) = | (3 log(z/Mf) (:c/Mf - 1) 1 + log (x/Mf,) ( x/M 2 H - 1) 1 ) 


(48) 


Mr and Mz denoting the Higgs and the Z boson mass, respectively. Hence, the overall 
light neutrino mass matrix, m v , is given by the sum of the tree-level ( [TT| and one-loop (47) 
terms, which in the basis of charged lepton mass eigenstates satisfies the relation 


m v = m/ ree + ml loop = Up MNS diag(mi, m 2 , m 3 ) Uj 


PMNS ' 


(49) 


The finite radiative correction given in (47) is in general subdominant in the case of RH 


neutrinos with a high mass scale M S> v, but it may be sizable and comparable to the 
tree-level term in seesaw scenarios where the lepton number violating scale is taken below 
the TeV range. It is therefore interesting to analyse in greater detail the dependence of the 


light neutrino masses on the additional finite one-loop contribution, Eq. (47). 
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In the basis in which the RH neutrino mass is diagonal, the one-loop correction of interest 
has the following form: 




(OVYik Ml 


3 log (Ml/Ml) , log (Ml/M, 


+ 


Hi 


(50) 


(4 n V? -* V MU - 1 MU Ml - 1 

where we have used Eqs. i§ and (Jr|) . The contribution of the one-loop correction under 
discussion to the effective Majorana neutrino mass generated by the light Majorana 

neutrino exchange, as can be shown, is given by 


= (™y loop ): e . 


(51) 


The Scheme with Two RH Neutrinos 


In the phenomenologically interesting scheme with two RH neutrinos, for each non-zero 


eigenvalue rn k of Eq. (49), we have the exact relation 


0 = det[m fc l 3 x 3 + m D M R l (1 2X 2 - U{M R )) m R ] 

= m k det [m k 1 2X 2 + M^ 1 ( l 2x2 - U{M R )) m T D m D ] , 


(52) 


where the second equality follows form the Sylvester’s determinant theorem and we have 
introduced the function 7 


U{Mr) = 


(47m) 2 


F(M r M' r ) + M R F(M' r M r ) m r ') . 


(53) 


Using (52) and (20), we get the identity 


det 11 


2x2 


H(M r )] mf ee mf e = mi m h , 


(54) 


where mi (mh) is the smaller (larger) non-zero active neutrino mass, whose experimental 
value in the cases of NH and IH neutrino mass spectrum is given in Eqs. ( Jl7| and (18), 
respectively. 8 Therefore, the determinant on the left hand side of Eq. ( ]54| provides a mea¬ 
surement of the deviation of the tree-level mass eigenvalues from the observed neutrino 
masses. Notice that, this is a positive quantity smaller than one in the scenarios considered 


7 The definition given in Eq. (53) is by construction basis independent. 

8 In the convention we are using ?n) ree m) t ree = m^ ee m^ ee ( mf ree m^ ree = and mimh = mim-s 

(rnimh = TO 1 TO 2 ) in the NH (IH) case. 
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here. As a consequence of Eq. (54), one has that in the case mf ee = 0, i.e. if two of the 


active neutrinos are massless at tree-level, it is not possible to generate at one-loop level two 
non-zero light (active) neutrino masses in the spectrum. In other words, in such a scenario 
both the solar and atmospheric neutrino oscillation mass differences cannot be radiatively 
generated. 

As it is not difficult to show, in the minimal scenario with only two heavy Majorana 


neutrinos, in which condition (28) is exactly fulfilled, the one-loop contribution to the light 


neutrino mass matrix goes to zero in the limit AM = M 2 — Mi —0. Indeed, from Eqs. (28) 


and (51) we find: 


m 


1—loop 

w 


(4 7 tv) 2 
— z(2 + z) 


{0V)h Ml 


3 log (Ml/M 2 Z ) , log {Ml/Ml) 


Mf /Mf 


+ 


Mf/M 2 - 


{Ml -)■ Mf) 


3 log(Mf/Mf) + log(Mf/Mf) 


Mf/Mf - 1 Mf /Mjj - 1 


(55) 


where z = AM/M 1; i.e., M 2 = (1 + z)Mi. Note that Eq. (55) is valid for arbitrary values of 
z and Mi. In the case of Mf, Mf <C Mf, Mjf we get: 

m 2 t 


m 


1—loop 

w 


lei 


(4 7Tu) 


Mf [8(1 + z) 2 log(l + z) + z(2 + z) (3 log(Mf/M|) + log(Mf/Mj))] 


(56) 


If, in addition, z -C 1, this expression further simplifies to: 
{6Vf 


m 


1—loop 


/el 


Mf z(2 + z) [4(1 + z) 2 + 3 log(Mf/Mf) + log(Mf/M|)] . (57) 


4/3 (4 7 m ) 2 

In the opposite limit, namely, Mf, Mf Mf, Mf), loop takes also a rather simple form 
for 1. In this case, to leading order in z <C 1, we obtain: 

1 


1 —loop o 

m 3fi = ~2z 


(Attv ) 2 


(3M 2 +M 2 ) . 


(58) 


Thus, in the scheme considered here, in which condition (28) is fulhlled, the magnitude 


of the one-loop correction to of interest, exhibits a strong dependence on 

z. This dependence is particularly important in the case when the two heavy Majorana 
neutrinos form a pseudo-Dirac pair, 0 < AM Mi,M 2 , or z -C 1. In this case the ratio 
of the one-loop correction to the 0 u/3/3 decay amplitude and the heavy Majorana neutrino 


1—loop 


exchange contribution given in Eq. (32), mf ; g loop /m*))) avy |, practically depends only on the 
mass Mi. As it is not difficult to show, for f(A ) = 0.79 (0.033), i.e., for 76 Ge ( 48 Ca), we have 
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1 10 10 2 10 3 10 4 
Mi(GeV) 


FIG. 2: The contributions to the Ou/3/3 decay effective Majorana mass due to the one-loop cor¬ 
rection to the light neutrino mass matrix (dashed line) and due to the heavy Majorana neutrino 


exchange (solid line), \fnpp°° V \ and \rn^ e g Vy \ (Eqs. (55) and (32)), as functions of the heavy Majo- 




rana neutrino mass Mi, for AM = 10~ 2 GeV, \(0V) e i\ 2 = 10~ 3 and f(A) = 0.079 (i.e., for 76 Ge). 
The range of values the effective Majorana neutrino mass can take in the case of light Majorana 
neutrino exchange and IH spectrum is also shown (the band in red color). See the text for further 
details. 

|m^ loop /mJ e ; vy | « 1 at M 1 « 15 (9.7) GeV. For M, > 15 (9.7) GeV (M, < 15 (9.7) GeV), 


W 

|m^ loop | is bigger (smaller) than |m^ avy 


. This is illustrated in Fig. 


which shows the 


dependence of |m^g loop | and m!)^ lvy on Mj > 0.5 GeV for AM = 10 2 GeV in the scheme 




in which condition (28) is exactly fulfilled and fixing the active-sterile mixing to the reference 
value of \(9V) 2 el \ = 10^ 3 . In this plot the Higgs mass has been set to M H = 125 GeV. Note, 
however, that given the values of M z = 90 GeV and M H = 125 GeV, for Mi = 15 (9.7) GeV, 
the factor (4(1 + z ) 2 + 3log(M 2 /M|) + log(Mf/M^)) in Eq. (57) for m^ loop is negative. 
Thus, at Mi = 15 (9.7) GeV, we have m i“ loop / m ^ avy > q (see Eq. (32)), and therefore a 


cancellation, or even a partial compensation, between the two terms m^ 10015 and m|) e g avv i 


"A9 


in 


the 0 u/3/3 decay amplitude is impossible. 

As it should be clear from Fig. [ 2 ] and Eqs. (56 58), |m^ loop | grows rapidly with the 
increase of M\. However, the dependence of |m^ loop | on z when z « 1 makes it possible, 
in principle, for |m^ loop | to have values in the range of sensitivity of the current and next 
generation of 0 u/3/3 decay experiments, i.e., to have |m^ loop | ~ (0.01 — 0.10) eV even for, 
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e.g., M\ = 10 3 GeV and the maximal value of |(0V)g 1 | = 10 -3 allowed by the current data. 
This requires, however, exceedingly small values of z, which lead to a subleading heavy 
neutrino contribution. Indeed, using the quoted values of Mi and \{0V) 2 el \, and taking 


into account that v = 246 GeV, it is not difficult to find from Eq. (58) that we can have 
0.01 (0.10) eV for z « 6 x 1CT 10 (6 x 10~ 9 ). Such a small value of z suggests 


\m 


1—loop I 

W 


a severe fine-tuning, but it can also be understood in the context of the ISS scenario as a 
technically naturally small value of the lepton number violating parameters of this model. 


In the analyses which follow we will not assume that Eq. (28) relating ( 9V) e i and (0V ) e2 
is satisfied. We will use only the phenomenological constraint on (6V) e 1 and (9V) e2 [TDHT7] . 


Notice, however, that for values of the Casas-Ibarra parameter |q| > 6 (see Eqs. (16), (24) 


and (25)), the relation given in Eq. (28) is effectively satisfied. 


B. One-loop Generalisation of the Casas-Ibarra Parametrization 


In order to make sure that we generate the correct light neutrino mixing pattern, it 
is useful to generalise the Casas-Ibarra parametrization introduced in the previous section 
including the one-loop correction to the neutrino mass matrix. Taking into account the 


expression (50) for (ml loop )^' in the basis in which the RH neutrino mass is diagonal, 


Eq. (49) takes the explicit form: 


(m v )u> = ~ (m D V ) 


£k 


M ; 1 


(4 7T V ) 


Mk ,' 3 log (Mj/Mj) + log {MU Ml, 


Ml/Ml - 1 Ml/Ml - 1 


= -(m n V)e k A k 1 (y T m I D ) te = ((7p MNS diag(m 1 ,m 2 ,m 3 ) (7f MNS )«.. 


(V T m T ' 


d)w 


(59) 


Hence, in analogy to the tree-level contribution, we have now 

(±z m~ 1/2 V^ MNS 6 V M A" 1/2 ) m" 1/2 V^ MNS 9 V M A" 1/2 ) T = RR T = 1. (60) 

Thus, we get the following expression for the heavy Majorana neutrino couplings in the weak 
charged current, or equivalently, for the active-sterile neutrino mixing, at one-loop order: 


ev = =F i U pmns rn 1/2 R A 1/2 M 1 


(61) 


In the numerical analysis reported in section IV we will make use of this parametrization of 


9V, with R given in (21) and (22), in order to include the one-loop corrections to the light 
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Mj(GeV) 


Mj(GeV) 


FIG. 3: Maximum value of the contribution to the Ozz/3/3 decay effective Majorana mass due to the 
heavy Majorana neutrino exchange |m^ avy | (solid thick line) for 76 Ge and AM = 10~ 2 GeV in the 
IH (left panel) and NH (right panel) case, including the following constraints: |m^ avy | A 0.5 eV 
and \{6V) e i\ 2 + |(#V) e 2 | 2 < 2 x 10~ 3 . The corresponding values of the contributions to the Oz//3/3 
decay effective Majorana mass due to the tree-level (dashed line) and one-loop correction (dotted 
line) to the light neutrino mass matrix, |m^ e | and |m^ loop |, are also shown. The range of values 
the effective Majorana mass can take in the case of light Majorana neutrino exchange and IH (NH) 
spectrum is shown in the red (blue) band. See the text for further details. 


neutrino masses and at the same time ensure that all the neutrino mixing parameters match 
with their experimental values. 

In Fig. [3] we illustrate the interplay between the contributions to the 0 u/3/3 decay ef¬ 
fective Majorana neutrino mass due to the heavy Majorana neutrino exchange, mj)^ avy |, 
the tree-level light neutrino masses, = |(m* ree )*J, and the one-loop correction to the 

light neutrino mass matrix, |m^ loo, ’| = |(m 3_loop )* e |, using the generalised Casas-Ibarra 
parametrization derived above. In particular, we have maximised |m^ avy | over the free pa¬ 
rameters of the model ($ 45 , 7 and the Dirac and Majorana phases of the PMNS matrix), in 
order to show the maximum heavy neutrino contribution to the process (solid thick line) as 
a function of Mi for AM = 1CT 2 GeV and fixing the already measured PMNS parameters 
and neutrino squared mass differences to the best fit values given in [26j. The Higgs mass 
has been set to Mh = 125 GeV. I 11 the plot we show the corresponding value of the separate 
contributions associated to the tree-level (dashed line) and one-loop correction (dotted line) 
to the light neutrino mass matrix. We also impose the following constraints: |m^ av> | < 0.5 
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eV and \(9V) el \ 2 + \(9V) e2 \ 2 < 2 x 10“ 3 . 

From Fig. [3] we conclude that for Mi < 1 GeV the one-loop correction is subleading 
for AM = 1CT 2 GeV, being the tree-level contribution the one responsible for the light 
neutrino mass generation. At the same time, in that region the heavy neutrino contribution 
to the 0 u/3/3 decay effective Majorana neutrino mass can be sizable and larger than the one 
from light neutrino exchange. According to the estimate given in Fig. [lj for M\ < 1 GeV 
there is no need of any enhancement of the active-sterile mixing with respect to the naive 
seesaw scaling in order to obtain a sizable ||. However, around Mi ~ 2 GeV, the one- 
loop correction starts to be of the same size as the value of the light neutrino contribution 
dictated by neutrino oscillation data. Indeed, this correction increases with Mi in such a 
way that in order to stabilise the light neutrino mass and mixing, a fine-tuned cancellation 
between the tree-level and one-loop correction is required. This is reflected in the fact that 
for Mi > 5 GeV the dotted and dashed lines merge. Therefore, as it is shown in Fig. [3j 
for 5 GeV < M\ < 1 TeV a sizable |m^ avy | can in principle be realised, but a fine-tuned 
cancellation between the tree-level and one-loop contributions to the light neutrino masses 
is also necessary. 

Note that the bound |m^ av> | < 0.5 eV imposed by us can be saturated for Mi < 100 
GeV. At Mi = 10 GeV, for instance, we have |m^ avy | = 0.5 eV for |(6M) e i| 2 + |(6M) e2 | 2 — 
0.8 x 10~ 4 , where we have used f(A ) = 0.079 corresponding to 76 Ge. For Mi > 100 GeV the 
maximum value of |m^ av ' | decreases with M x since an active-sterile mixing \{9V) e i\ 2 bigger 
than 2 x 10~ 3 would be required in order to saturate the bound. 

It is interesting that the solid line and the blue and red bands in Fig. [3] intersect around 
Mi ~ 10 3 GeV. This implies that in the case of NH neutrino mass spectrum, the effective 
Majorana neutrino mass \mpp | can be larger at 0.1 GeV < Mi < 10 3 GeV than that 
predicted in the case of the light neutrino exchange mechanism. In particular, it can be 
in the range of sensitivity of the experiments aiming to probe the range of values of the 
effective Majorana mass corresponding to the IH and quasi-degenerate (QD) light neutrino 
mass spectra (see, e.g., [36|). In the case of the IH light neutrino mass spectrum, the 
indicated result implies that at Mi < 10 3 GeV there can be, in principle, a significant 
interplay between the light and heavy Majorana neutrino exchange contributions in the 
effective Majorana mass, as discussed in detail in [7] and summarised by us at the end of 
subsection II.A (see the paragraph before the last in subsection II.A). More specifically, due 
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to this interplay of the light and heavy Majorana neutrino contributions, | mpp\ can be larger 
(smaller) than that predicted in the case of the exchange of light neutrinos with 1H mass 
spectrum and \mpp | will exhibit a dependence on the atomic number A of the decaying 
nucleus. It should be mentioned that, given the already high level of fine-tuning required for 
the cancellation between the tree-level and one-loop light neutrino contributions in mpp, an 
additional cancellation between the light and heavy Majorana neutrino contributions would 
suggest further fine-tuning. 

The main features of Fig. [3] also appear for larger splittings AM. In particular, the 
necessity of fine-tuned cancellation between the tree-level and one-loop correction to the 
light neutrino mass matrix is present also in this case. The level of the fine-tuning required 


increases with Mi, as we will show in Section IV 


C. Radiative Corrections to the ESS and ISS Scenarios 


In this section, we compute the one-loop contribution to the effective Majorana neutrino 
mass in the ESS and ISS limits of the seesaw Lagrangian (JTj) with two RH neutrinos. Ac¬ 
cordingly, we apply the parametrization of the Dirac and Majorana mass matrices reported 


in Eq. (33) to the general expression given in Eq. (47). The exact result of the one-loop 


contribution in terms of the parameters introduced in (33) is reported in Appendix |A] 
For the ESS scenario we have at leading order in A/// 


m 


1—loop 

PP 


I 1 ' vie 


2 (4 


7T 


3 In (/r /2 /M|) In (// 2 /M 2 ) 




+ 


^/Ml 


Notice that for // 3> M H , M z , this expression reduces to 


m 


1—loop 


Vie 




(62) 


(63) 


m (4 jt) 2 

Therefore, when // Mh, M z , since the lepton number violating scale // is introduced at 
high energies, the one-loop contribution to the light neutrino masses appears to be suppressed 
as 1///, as expected. 

In the ISS realisation, i.e. for ev, /i, // <C A, we obtain 


m 


1—loop 

PP 


( * /i 9 \ /3 In (A 2 /My) In (A /Ml) 

(e A 2/ie 2 / 2c - 2 Vie) y \ 2/ M 2 _ + A 2/ M ^ _ 


(4tt) 


(64) 


li + li’ y\ e ( AMffMz — A 2 (Mfj + 3M|) In (A 2 /M 2 H ) 


(4 vr) 


(A 2 - Mf) (A 2 


MJ r ) 


+ 


3 In (A 2 /Ml) 
(A 2 /Ml - l) 2 ' (A 2 /Mf - l) 2 


+ 
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It is remarkable that in the ESS limit with // < Mu, Mz and in the ISS limit the one- 
loop correction to the light neutrino masses has a contribution proportional to /T. This 
dependence on fx' is very relevant since at one-loop the light neutrino contribution to the 


0 v(3(3 decay amplitude does depend directly on //, as for the heavy contribution in (41) 


and (45). This makes much more difficult to obtain a dominant contribution from the 


RH neutrinos in this limit, unless a fine-tuning of the seesaw parameters is introduced to 
guarantee the smallness of the neutrino masses as it was indeed already shown in Fig. [3j 

IV. LARGE HEAVY NEUTRINO CONTRIBUTION TO 0v(3/3 DECAY 

In this section, we will address in more detail the question if the RH neutrinos can 
eventually give a sizable contribution to the 0zz/3/3 decay rate. As we have already mentioned, 
cosmological constraints close the mass window of M < 100 MeV and thus only if the 
RH neutrino masses are larger than 100 MeV, a direct contribution to the process of interest 
can be expected. 

Following the notation in Ref. [2], the 0 v/3(3 decay rate can be written as 

2 


L 0 v0/3 

In 2 


= G, 


01 


/ / bp 


(65) 


where Goi is a well-known kinematic factor, U is the unitary matrix given in Eq. Q which 
diagonalizes the complete neutrino mass matrix, rrij are the corresponding eigenvalues, i.e., 
the neutrino masses (light and heavy), and Af 01 ^ are the Nuclear Matrix Elements (NMEs) 
associated with the process. Notice that the NMEs depend on the mass of the neutrino 
mediating the process since the dependence on the neutrino propagator is already included 
in the NMEs computation. The sum should be made over all the neutrino masses, including 
the heavy ones. In the following we will use the NMEs data provided in [2]. In particular, 
we will consider the NMEs computed for the ' 6 Ge. However, we have checked that the 
conclusions of our analysis do not significantly change considering a different nucleus. 

We will use the modified Casas-Ibarra. parametrization of the active-sterile neutrino mix¬ 


ing given in Eq. (61), to compute the full effective Majorana neutrino mass mpp, which is 
given by the sum of the contributions from the exchange of the light and heavy Majorana 
neutrinos. In this way, we include in the computation the effect of the one-loop correction 
to the light neutrino masses, reproducing at the same time the correct neutrino oscillation 
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parameters. We will also take into account the relevant bounds on the active-sterile mixing 
which come from direct searches, charged lepton flavour violation and non-unitarity con¬ 
straints HDHIZ!. Notice that the inclusion of such bounds guarantees the perturbativity 
of the neutrino Yukawa couplings for any value of RH neutrino masses considered in this 
paper. 

In the top panels (down panels) of Figs. [I] and [5j the blue shaded area corresponds to the 
region of the parameter space in which 1CT 2 eV < + m^g av> | < 0.5 eV (10 -2 eV < 

I < 0-5 eV), projected on the 7 — AM plane for NH (IH) and several values of M\. In 
these plots we have fixed the already measured PMNS parameters and neutrino oscillation 
mass differences to the best fit values given in [26]. The relevant Majorana and Dirac CP 
violation phases in the PMNS matrix have been set to zero, but we have checked that 
there is no significant impact on the results when other values are considered. The Casas- 
Ibarra parameter d 45 is also set to zero. It is irrelevant when the heavy Majorana neutrino 
exchange contribution is dominant (subdominant) in rripp , but can play an important role in 
the interplay of the light and heavy Majorana neutrino exchange contributions when these 
two contributions are comparable in size [?]. The Higgs mass has been fixed to Mh = 125 
GeV. The solid black line stands for different values, stated in the plots, of the a parameter 
defined as 

a = KMl/Kfl. (66) 


where f = m^ e + rri^ loop is the full (tree-level plus one-loop) contribution to 171 / 3/3 given 
by the light neutrinos. Therefore, a quantifies the level of fine-tuning in the cancellation 
between m jj r J e and rn^ loop described in section III B and required in order to keep the light 
neutrino masses and mixing to the observational values. Notice that the level of fine-tuning 
increases with a. The region to the right of the black solid line corresponds to values of a 
larger than those stated in the plots. 

In the red shaded area of Figs. [4] and [5j the ratio between the leading and next to leading 
order contributions to the light neutrino masses in the seesaw expansion is smaller than 10 . 
The next to leading order contribution is given by 


8m v = -- (mf 6 + m l ° op ) {6 V){0V)i - ^ (0V)* { 6 V) T (m‘ ree + m l ° op ) . (67) 


From this expression, one can conclude that a cancellation between the one-loop and tree- 
level contributions to the light neutrino masses remains at next to leading order in the seesaw 
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FIG. 4: Neutrinoless double beta decay (Mi > 100 GeV). The blue shaded areas 

in the top panels (down panels) represent the region of the parameter space in which we have 
10“ 2 eV < |mjg ht + mj^ avy | < 0.5 eV (10“ 2 eV < |mj^ avy | < 0.5 eV) in the case of NH (IH) 
neutrino mass spectrum with the active-sterile mixing (or couplings) (OV)ik satisfying the bounds 
form direct searches, charged lepton flavour violation and non-unitarity constraints. The black 
solid line stands for different values of the parameter a = | m i~ loop |/| m ^g ht | ; w hi c h quantifies 
the fine-tuning required in order to achieve the cancellation between the one-loop and tree-level 
contributions to the light neutrino masses. In the region to the right of the red dashed line the 
ratio between the leading order and the next to leading order contributions to the light neutrino 
masses in the seesaw expansion is smaller than 10. The gray region to the right of the dotted line 
corresponds to yi e r n 1 pp 00 ' c ' > 16 7r 2 m^ ht . The blue dashed line corresponds to |m^ avy | = 0.05 eV. 
The measured neutrino oscillation parameters are fixed to the central values reported in [26]. 


expansion. This is in agreement with Figs. [4] and [5j which show that the next to leading 
order contribution is always negligible in the range of parameters of interest. 

Ignoring for the time being the impact of the two-loop corrections, which will be com- 
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merited below, two main conclusions can be extracted from Figs. [| and [5j First, we have 
proved that a sizable and dominant heavy neutrino contribution to the 0 u/3/3 decay is pos¬ 
sible for RH neutrino masses as heavy as 10 TeV, satisfying at the same time the relevant 
constraints and keeping under control the light neutrino mass and mixing pattern. Second, 
and not less important, it is shown that this possibility can only take place if a highly 
fine-tuned cancellation between the tree-level and one-loop light neutrino masses is at work. 
The level of fine-tuning ranges from a = 10 4 to 10 9 , for heavy masses between Mi = 100 
GeV and Mi = 10 TeV. On the other hand, the level of fine-tuning is smaller for lighter 
masses, being in the case of Mi = 100 MeV smaller than a = 2. In addition, we have 
checked that for Mi > 10 TeV a heavy contribution to mpp in the range of sensitivity of the 
next-generation of experiments, \mpp\ > 0.01 eV, cannot be expected. 

Figs. [4] and [5] also show that in the limit AM Mi the sizable heavy neutrino contri¬ 
bution corresponding to the blue region becomes independent of AM, according with the 
ESS limit - see Eq. (32). However, in the ISS limit AM <C Mi this is not the case and, 
according to Eq. (32), the smaller the heavy splitting AM, the larger is the value of 7 . 

Notice that in the IH case we have plotted only ni^T' y because | is already in the 


planned range of sensitivity of the next generation of 0 u/3/3 decay experiments. In this case 
for Mi < 10 3 GeV and AM << M ])2 , there can be, in principle, a significant interplay 
between the light and heavy Majorana neutrino exchange contributions in the effective 
Majorana mass, as discussed at tree level in detail in [7] and summarised by us at the end 
of subsection II. A (see the paragraph before the last in subsection IF A). More specifically, 
due to this interplay of the light and heavy Majorana neutrino contributions, \mpp\ can be 
larger (smaller) than that predicted in the case of the exchange of light neutrinos with IH 
mass spectrum and \mpp | will exhibit a dependence on the atomic number A of the decaying 
nucleus. This can happen roughly in the region located to the left of the blue dashed line 
corresponding to |m^ av ' | = 0.05 eV inside the blue areas in Figs. |4]and 5 

In the NH case, the light neutrino contribution is smaller than 10 ~ 2 eV and therefore any 
sizable effect to the process is due to the heavy neutrinos. This is why in the NH case we 
plot the total contribution mpp, including light and heavy neutrinos. 

It follows from Fig. 4 that for Mi > 100 GeV the regions of interest (the blue shaded 
areas) correspond to 7 > 6 . For such values of 7, as it is not difficult to show, we have for 
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FIG. 5: Neutrinoless double beta decay (Mi < 100 GeV). The same conventions as in 
Fig-B but for different choices of M\. 


the NH and IH neutrino mass spectra: 


M 2 \(6V) e2 \ 


Mi \(9V) el \ 

d 2 7 

T 

; 2 7 

T 

Taking into account that Fig. 4 is obtained by setting to zero the phase 9 45 and the Dirac 
and Majorana phases in the PMNS matrix and by using the best fit values of the neutrino 


g ' 2 

~ —7- | U e2y /m^~ iU e3y /m^\ 


~ \U elx /m^~ iU e2 Jm 2 \ 


NH 


IH 


( 68 ) 

(69) 


oscillation parameters, Eqs. (68) and (69) imply the following relations between | (6 > V) e i( e2 )| 2 
and the parameter 7: 


M^ev) 


el 


M 2 \(6V) e2 1 2 w e 27 0.94 (12.4) x 10“ 3 eV, NH (IH). 


(70) 


In view of the high level of fine-tuning required in order to have a cancellation between 
the tree-level and one-loop light neutrino masses, the obvious question arising here is what 
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is the role of the two-loop corrections. Can the two-loop corrections spoil this fine-tuned 
cancellation? In order to answer this question, we estimate the impact of the two-loop 
contributions. Since we are studying the case in which heavy neutrinos can give a sizable 
contribution to the Ou/3/3 decay, which means relatively large Yukawa couplings, we expect 
the diagram with two Higgs bosons in the loop to be the leading two-loop contribution to 
the light neutrino mass matrix. The contribution of this diagram can be roughly estimated 
as 


7 .2 

2—loop y le 1—loop 

m w (i 


( 71 ) 


where m^ loop is the one-loop contribution in This estimate of the impact of the 

two loop corrections is also shown in Figs. |4] and |5j where the gray area to the right of the 
dotted line corresponds to the region of the parameter space with y\ e m^ loop > 167r 2 m^® ht . 
This region of the parameter space is excluded since the two-loop correction, which would 
dominate the light neutrino masses, would be larger than the value dictated by neutrino 
oscillation data. Notice that this would essentially exclude the possibility of having a large 
sterile neutrino contribution for M\ > 1 TeV, as can be seen in Figs. [4] and [5} For Mi < 100 
GeV the impact of the two-loop correction is basically negligible. 


V. CONCLUSIONS 

We have performed a systematic analysis of the radiative corrections to the light neutrino 
masses arising in low scale type I seesaw scenarios, where the RH (sterile) neutrino masses 
vary in the interval 100 MeV < M < 10 TeV. Within this range of masses a significant 
enhancement of the neutrinoless double beta (Ou/3/3) decay rate in several isotopes - at the 
level of sensitivity of the present and next generation experiments searching for this rare 
process - is possible, due to the new physics contribution in the decay amplitude given by the 
exchange of the virtual heavy sterile neutrinos. Notice that one of the most clear signatures 
of a significant heavy sterile Majorana neutrino contribution to the 0 v(3(3 decay amplitude 
is the dependence of the 0 u/3/3 decay effective Majorana mass, \mpp |, on the atomic number 
A of the decaying nucleus jUj. 

The requirement of a sizable contribution of heavy neutrinos with masses > 1 GeV to 
the 0 v/3/3 decay implies strong cancellations between the tree-level and one-loop expressions 
in the light neutrino mass matrix m v originated from the seesaw mechanism. We show that 
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such a cancellation can always be achieved while being consistent with neutrino oscillation 
data and low energy constraints from direct searches, charged lepton flavour violation and 
non-unitarity by using a generalisation of the Casas-Ibarra parametrization of the neutrino 


Yukawa matrix, which can be derived from Eqs. (60) and (61). We clarify the connection be¬ 
tween this parametrization and the lepton number breaking terms in the seesaw Lagrangian, 
as usually defined in extended as well as inverse/direct seesaw UV completions of the Stan¬ 
dard Model. Then, we numerically quantify the level of fine-tuning between the tree-level 
and one-loop parts of m v in the case the heavy neutrino contribution m^ avy to the effective 
Majorana neutrino mass - which enters in the 0 u/3/3 decay amplitude - is sizable, namely 
K e ; vy | > 0.01 eV. 

The main results of our analysis are summarised in Figs. [4]and[5j where we show that a 
fine-tuning of one part in 10 4 (10 5 ) for RH neutrino masses ~ 100 (1000) GeV is unavoidable 
in order to have an observable effect in 0 v(3/3 decay experiments. Furthermore, we conclude 
that for seesaw scales M larger than few TeV, two-loop effects in the generation of the 
light neutrino masses cannot be neglected, thus excluding the possibility of having a large 
|m^ avy [. Conversely, in the low mass regime, M < 1 GeV, the level of fine-tuning in the 
seesaw parameter space is very mild and the sterile neutrino contribution can easily exceed 
the current limits on the effective Majorana neutrino mass. 

Finally, we can conclude on the basis of the results obtained in the present analysis 
that 0 v(3f3 sets the strongest constraints on lepton number violation in low scale type I 
seesaw extensions of the Standard Model. In particular, this implies a strong suppression 
of processes which involve the production at colliders (LHC included) of RH neutrinos and 
their decays with two like-sign charged leptons in the final state (see, e.g., M). 
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Appendix A: m* loop in an arbitrary basis for 2 RH neutrinos 


We report in this appendix the full computation of the one-loop correction (47) to the 


light neutrino mass matrix in terms of the seesaw parameters introduced in Eq. (33), from 


which it is possible to derive the one-loop correction to effective Majorana neutrino mass in 


the extended and inverse seesaw limits, Eqs. (62) and (64), respectively. In order to obtain 


an analytic expression for the one-loop neutrino mass matrix, we conveniently change the 
basis of the heavy RH neutrinos, i.e. u aR = V ab is' bR , with the unitary transformation 

(Al) 



In the new basis the RH neutrino Majorana mass matrix takes the form: 

. ~ rp ~ 1 / 2A — (/! + //) —*(// — //') 

m' r = v t m r v=-\ 

z y -i(n-n') 2A + (fj, + //). 

Then, the resulting one-loop Majorana mass term for active neutrinos is 


(A2) 


m 


1—loop _ 


( 47 ™)' 


m D V ( M' r ' F(M' r M'J) + F(M' r 'M' r ) M' r ‘ ) V T m T D , (A3) 


r 't 


7 t 


ff 


where the loop function F(x) is defined in Eq. (48) and the Dirac mass matrix me is 
parametrized as in (33). In this case we have: 9 


M' r M' r — ^A^ + - (|/x | 2 + Ip'l 2 )^ ( 12 x 2 — A(a, b, c)) , 

M'r Mr = ^A 2 + - (|p | 2 + |p 7 | 2 )^ ( 1 2X 2 - A(a,- 6 , c)) , 


where 


(A4) 

(AS) 

(A6) 


A(a,b,c) = (acr 3 + bcr 2 + ccri) , 

cTj (z = 1, 2, 3) denoting the 2x2 Pauli matrices. The real parameters a, b and c are defined as 

2 A Re (n + fi') 


a = 


b = 


c = 


2 A 2 + H 2 + |/i '| 2 

2 A 2 + |p | 2 + |p '| 2 
2 A Im (, n 7 — [/,) 


2 A 2 + 


+ \A 


(A7) 

(AS) 

(A9) 


We assume without loss of generality that the parameter A in (33) is real. 
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In this way, one can obtain a closed form for the logarithms which enter in Eq. (A3) through 
the loop function. Indeed, we have 


with 


log [l 2x 2 - A(a,b,c)] 


Y -A(a,b,c) 

n=l 


(A10) 


[A(a,b,c) n } u — t ( a 2 3 + b 2 + c 2 ) 1+,t| (a (1 - (-1)”) + Va 2 + bA + c 2 (1 + (-1)")) , 

[A(a, 6 ,c )"] 22 = -I (a 2 + b 2 + c 2 )='- 1+ ”> (a (1 - (- 1 )") - Vo 2 + V + c 2 (1 + (- 1 )")) , 
[A(a,M ”] 12 = p(a, &,<:)% = (6 + ic) (a 2 + 6 2 + (-1 + (-1)") . (All) 


Then, one can show that the infinite series in (A10) gives the exact results 


r , , n 2 a tanh 1 ( \/a 2 + b 2 + c 2 ) . „ „ ox 

log [l 2 x 2 - A(a, 6, c)] = - y -log (l — a -6 -c) , 

ya + 

[log[l 2X 2 -A(a,b,c)\] 12 = [log(l 2x2 - A{a,b,c)))\y = 

r , 2 a tanh -1 (\/a 2 + 6 2 + c 2 ) , 

[iog[i 2 x 2 -A(a,6,c)]] 22 =- 7a2 + y2 + g2 -i°g(i 




a /a 2 + b 2 + c 2 

i 2 -b 2 - c 2 ) . (A 12 ) 


Therefore, by replacing Eqs. (A 12 ) in (A3), we obtain an analytic expression for the one-loop 
contribution to the light Majorana neutrino mass matrix as a function of the parameters 


given in (33). 
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